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Dynamics of hadro-quarkonium system is formulated, based on the channel coupling of a light 
hadron (h) and heavy quarkonium (QQ) to intermediate open-flavor heavy- light mesons (Qq, Qq). 
The resulting effective interaction is defined by overlap integrals of meson wavefunctions and (hqq) 
coupling, where h is -n,p,uj,(f>, without fitting parameters. Equations for hadro-quarkonium ampli- 
tudes and resonance positions are written explicitly, and numerically calculated for the special case 
of nT(nS) (n = 1,2,3). It is also shown, that the recently observed by Belle two peaks Z(,(10610) 
and Zj,(10650) are in agreement with the proposed theory. It is demonstrated, that theory predicts 
peaks at the BB* , B* B* thresholds in all available 7rT(n5') channels. Analytic nature of these 
peaks is investigated, and shown to be due to a common multichannel resonance poles close to the 
BB*,B*B* thresholds. The general mechanism of these hadro-quarkonium resonances does not 
■ assume any molecular or four-quark (tetraquark) dynamics. 
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I I. INTRODUCTION 

&■ 

It was found in experiment [l|-[4| that resonances may appear in the system of a hadron and heavy quarkonium, 
which may be called hadro-quarkonium, see [f| for a review. On theoretical side the prevalent approaches associate 
hadro-quarkonia with molecular or four-quark (4g) states [6|-[16|. In the first case hadro-quarkonia are weakly bound 
states of two heavy-light mesons of the closest threshold with interaction tuned to produce loosely bound or virtual 
states, and in the 4<? states thresholds cannot be easily connected with 4c/. However, it will be argued that channel 
I/"") ' coupling (CC) near thresholds may play the dominant role in hadro-quarkonium dynamics, as was shown for heavy 
lO ', quarkonia in our previous papers [^Oil 1 . It was shown there, that strong CC, calculated basically without fitting 
parameters, shifts the 2 3 Pi (cc) pole exactly to the D D* threshold. In this way the A (3872) phenomenon was 
explained using only one parameter M u , which was fixed in previous studies [24j-|27| and universal input: the string 
tension a, the current (pole) quark masses, and the strong coupling a s (q). Recently M u was found from the first 
principles in QCD [28]. It was shown there, that M w can be calculated as the matrix element of the operator err, 
where a is the string tension and r is the length of the string. The decay width of "0(3770) is reproduced in this 
way and corresponds to M w ss 0.8 GeV. Our starting point is the first principle derivation of the CC interaction of 
standard heavy quarkonia with open flavor channels, using strong decay theory [28]. Similarly, in the phenomenon of 
A(3872), the systems wJ/ip and pJ/ip may take part with the thresholds near those of (DD* + h.c.) states. In the 
same way additional pions in the decay vertex appear with the only extra factor in the denominator f„ — 93 MeV. 
As will be shown below, the strong interaction of pions with Y(nS) mesons produces charged Z-type resonances. 
Recently in a series of papers 0, HH the CC methods have been successfully applied to the transitions in systems, 
containing heavy quarkonia and pions, or r\ meson and in this way the main features of experimental pionic spectra 
in reactions X' — > X tttt, were explained, together with kaonic and 77-meson final states (2fj|. |27|. Below we extend the 
formalism of channel coupling (CC) developed in [l7l [l8j to the case of a hadron h = n, (f>, 77, p, ui, ... interacting with 
the QQ state. 

We study the interaction and possible poles of hadro-quarkonium amplitudes in the formalism of (l7| . We assume, 
that the most important interaction in hadro-quarkonium is due to intermediate states of heavy- light mesons (Qq) (Qq) 
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1 A similar in spirit, but different technically the so-called rescattering model was developed and applied in particular to dipion transitions 
in bottomonia in [l9l423l ] 
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FIG. 1: The chain of transitions of hadro-quarkonia (h + QQ) and pair of heavy-light mesons (Qq) (Qq). h denotes light hadrons 

(h = TT,<p, ...). 



(e.g. DD, DD* , D* D* , D S D S , ... in case of hadro-charmonium). Therefore one should sum up the whole series of 
bubbles, consisting of hadro-quarkonia and heavy-light mesons, as shown in FigfT] To find poles in amplitudes, one 
can start and finish with any state, since poles belong to all channels, i.e. h (QQ) and (Qq) (Qq). We shall formally 
study the amplitudes for the transition of two heavy-light mesons again into the same or other pair of heavy-light 
mesons. 

It is important to stress, that in our mechanism of hadro-quarkonium resonances there is no direct interaction 
neither in the hadron-quarkonium channel, nor in the channel of two heavy-light mesons. The only interaction, 
which generates resonance poles, is the CC interaction, transforming hadro-quarkonium system into double heavy- 
light system. Therefore hadro-quarkonium resonances, predicted in our theory, is a clear example of CC resonances, 
introduced and calculated earlier in [2{|. We show below, that direct molecular resonances of BB* (if any) arc 
displaced and splitted in hadro-quarkonium in different hadro-quarkonium channels. 

To find the poles, we can use the so-called Weinberg Eigenvalue Method (WEM), discussed in detail in [l7|. It 



allows to define not only poles, but also resonance wave functions and was successfully applied in 17 1 to charmonia 

" l3 that X(3872) 



states, and in particular to X(3872), in situation of strongly coupled channels. It was shown in [Tj 
is due to bare n = 2 3 Pi resonance, shifted exactly to the D^D^ threshold by CC and the detailed experimental form 
was reproduced in [Tt], EH with a tiny cusp at D + D*_ threshold and no other bumps. No connection to pJ/ifi and 
cj J channels was taken into account in [171 H8| . assuming the corresponding partial widths to be generally small, 
and here we establish formalism for these channels, and 4> J /''Pi which allows to find out, whether the CC interaction in 
these cases is strong enough to produce poles. The situation with the uj J/ip channel is especially interesting, since it 
contains the resonance of its own, F(3940) [J, but in addition the decay X(3872) — > u> J/tp, found in [3(j, (see [3l| for 
a recent review) suggests that the whole CC system for JT(3872) should contain channels 2 3 Pi (cc ), DD* , lu J/ip and 
p Jhb. The CC analysis of this system in another framework (Resonance-Spectrum Method) (32|, was done recently 
in [33j. 

A special case of hadro-quarkonium is the pion-quarkonium system, where the CC interaction vertex is proportional 
to 1/ fn and numerically large, which might support the appearance of ir(QQ) resonances. Below we shall study 
specifically the case of tttt transitions in the T(nS) states, where these resonances appear in the final states T(n'S) irir. 

The plan of the paper is as follows. In section 2 some basic equations of WEM in the hadro-quarkonium case are 
written, and in section 3 those are exploited to write down exact equation for the possible poles in the general case of 
three sectors. In section 4 the special case of pion-quarkonium system is treated in detail and (nT(nS)) are found for 
n = 1, 2, 3. In section 5 results of calculations are given, and section 6 contains conclusions, comparison of molecular 
and CC dynamics and outlook. Two appendices are devoted to detailed derivation of decay transition kernel and the 
form of wave functions. 



II. DYNAMICS OF STRONG CHANNEL COUPLING FOR HADRO-QUARKONIUM 

We consider two strongly interacting sectors: sector I with heavy quarkonium state (QQ) plus hadron h = tt,lu, p, 77, <j) 
etc., and sector II, consisting of two heavy-light mesons (Qq) (Qq), in case of hadrocharmonium, it could be 
DD, DD*, D*D* , DD U D S D S etc. 

It is important, that we neglect interaction between any white objects, considering the limit of large N c . It means, 
that in our treatment there is no direct interaction between hadron and heavy quarkonium, as well as between heavy- 
light mesons. Justification of this approximation can be found in the fact known from NN interaction, that the main 
part of long-range forces between white objects comes from the exchange of one pion or a pair of correlated pions, 
which in case of deutron yields a small binding energy. However heavy quarks in heavy-light mesons do not contribute 
in this process and hence one-pion exchange in the system of two heavy-light mesons should be much smaller, that in 
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FIG. 2: The diagram of the hadron interaction V 



the NN system. This is also supported by the fact, that AiV, SAT and AA, EE interactions are relatively weaker, that 
the JVJV interaction. Therefore from our point of view in the molecular models of exotic charmonia one should take 
into account that much stronger attraction near the threshold occurs due to CC interaction between sectors I and II. 
One more support of this comes from our recent study of A (3872) dynamics in frR Il8j . where we have shown, that 
CC alone strongly shifts 2 3 Pi cc level by ~ 60 MeV to the D Dq threshold at 3872 MeV. 

As was shown in 17] . to study dynamics of CC in our case, one can reduce problem to the one-channel case, where 
another channel enters via the CC interaction V a ba, and a, b refer to sectors I, II respectively. If one is interested only 
in the possibility of bound states or resonances due to CC, one can start with any channel, and we shall work mostly 
in channel II and consider the amplitudes shown in FigfT]which are generated by the interaction ^212- This interaction 
in the formalism developed in [17| can be written in momentum space as the amplitude of the loop diagram, shown 
in Fig(2J with hadron (h) and quarkonium (QQ) in the n-th state 
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(1) 



with E n (jp) — \/p 2 + M 2 , where M n is the position of the n-th bare quarkonium (QQ) state. Indices 712^3, n' 2 n' 3 

and the overlap 



denote in- and out- quantum states of heavy-light masons, the hadron energy is Wft(k) = Wk 



•ir 



matrix elements Jnn 2 n 3 (p, k) define the probability amplitude for the transition of two heavy-light mesons (Qq) n2 , 

(Qq)n 3 with quantum numbers n2,«3 to quarkonium n-th state (QQ) n plus hadron h. One can derive Jnn 2 n 3 as a 
matrix clement of a hadron emission operator between wave functions of quarkonium and two heavy-light mesons 
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+ q)^7(q)C„ 3) (q-k) 



(2) 



where N c is the number of colours, r = c (u — v) 



We point that the w.f Vn" 2) , ^" 3) in © are no 



longer full w.f. of mesons, but the radial part R^qq-> ^Qq \ divided by V^r, while the angular part of the w.f. 

is accounted for in the factor 2/123. This transition kernel 2/^3 contains a coupling constant g q h of hadron with quark 
pair (qq), entering the hadron string-breaking (hqq) Lagrangian 



e lkx d A x 



(3) 



and another part, which comes from the Dirac trace of 7 matrices corresponding to the vertices in state (QQ) n and 
(Qq)n 2 , (Qq) 713 • 

To obtain the full vertex 1/123 m (2), one can use either the (4 x 4) form given in [24l HH or else the (2 x 2) form for 

wave functions and vertices, introduced in [13], Appendix B. Exact expressions for are given in Appendix 1 for 
the convenience of the reader. In this way for the n\ state of quarkonia and vector hadron one can write similarly to 
(B.3) 



(4) 
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and 

r<&{D) = -L, r% d (D*) = ^L, r^(i--(QQ)) = ^|. (5) 

Therefore in the case 1 (QQ) + h — > DD* , one obtains 

V123 = i 9h ej e ijk (6) 

where e is the polarization vector of a hadron and e^k is the Levi-Civita symbol. 

For Nambu-Goldstone bosons (-7T, K, n) the transition kernel was obtained in a different way in [24] . Indeed, pions 
accompany string breaking, yielding a coefficient 75 instead of g^e in ([3]), where M u is calculated via string tension 
a in [28], M u — 0.8 GeV. This is used in section 4 below, details are given in appendices 1 and 2. 

To define possible resonance position and wave function it is convenient to use the Weinberg Eigenvalue Method, 
which was extended to the case of coupled channel problem in [l7j . The corresponding equation for eigenfunction 
^f v (r,E) and eigenvalue rj„(E) can be written as 

H 9 v (t,E)+ f V ^ V ^ E ) ^(r',£)rfV = E* v (r,E) (7) 
J ilu{E) 

with boundary condition ^ u (r — > 00, E) ~ exp(ikr)/r, ^f v (0,E) = const and index v labels the discrete eigenvalues 
and eigenvectors. In the momentum space one can write 

vw 3 (p,£) - -^4y E / p G l(p)4,» s „;(p.p'«)i>i(p') 

At this point one realizes^that Eq.® can be seriously simplified, using the structure of the overlap integrals in ([2]). 
Indeed, it was shown in [24], H^] that wave functions of heavy-light mesons D,D* (B,B*) can be represented by the 
Gaussian functions with accuracy of the order of few percent. In this case the integral in ^ factorizes (see Eq.(2I) 
in[13) 

41 3 (P. k) = ( k ) xlt 3 (P) • (9) 

Moreover, it appears, that y4^„ 3 (k) are almost identical for the first two states of heavy-light mesons (e.g. B, B*) 

and are very close for the next two states (e.g. B Sl B*), hence we put by simplify (p^ n3 to ip( h \ 
As one can see in (fTJ), the integral on the r.h.s. also factorizes, when (0) is used, and one can write 

v£i inW3 (p,p', E ) = -E w x ^ ip) ^^ p ' )KM ' (10) 

n c 

where notation is used 

K(E)= [ d3k ^(k)^(-k) 
Insertion of (|I0l) in ([5} immediately yields 

^ns(P^) = W c Gi °tjE) E) E^^)^Us(P) E / |^^^(P')^^(P'^) (12) 

™ "2 "3 

which finally leads to a system of algebraic equations 



A n (£) = E/ pxit 3 (p)^« 3 (P,E), 



5 



where we have defined in (fT3 



So we obtain 



(nn'(E) — ^ y ^n2n 3 (P^) Xnn 2 n 3 (P) Xn'l l2n3 (p) ■ 



(14) 



det 



r h ,(E)6 nn ,-K n (E)C nnl (E) =0 (15) 



and one can look for poles E p , setting r\ v (E = E p ) = 1. 

Equivalently one can define amplitude in the sector I, h + (QQ), in which case the interaction "potential" V121 
assumes the form 



n 2 n 3 v 

and the WEM equation in sector I is 



d 3 p Jnn 2 n a (p, k) J^ 2 „ 3 (p, k') 



jB _ jff ^"3)( p ) 



(16) 



(k ' s) = - / OTO ^ ( k ' *> (k ' ^ (17) 

One can see from (fTB"]) , that is attractive and real for E below the DD (DD* ) threshold; and Gffl is 

Gi ° HKE) = H^TE = E n M+l h (k)-E- (18) 
The total Green's function in sector I has the form (see [ljj for discussion and details) 

G <-Hi^) = E ^V'1*| (2 ' £) < 19 > 

and near the resonance rj v (E) has the form 

Tfc,(£) = 1 + 77' (So - £ ) f S - E + ^) + ... (20) 



2 ' V 2 , 

Finally, one can define the i-matrix 

^y-yGy (21) 

which in the WEM can be written as 
where 

o„(k,£0 = (^o(k) - £)*W(k,£) (23) 



III. A GENERAL CASE OF THREE COUPLED SECTORS 



Till now only the connection of given channel h + (QQ) n to (Qq) (Qq) was considered. A more interesting situation 
can occur, when one adds also excited channels of (QQ) n > ■ An example of the physical situation of this kind is given 
by the 2 3 Pi(cc) state connected to ui J/i/j via the DD* channel. Hence in this case one has to consider three sectors: 
as before sectors I and II refer to h (QQ) n and (Qq) (Qq) respectively and sector III refers to (QQ) n ' states. One 
again writes equation for the wave function in sector II as in Eq. ([8]). but the interaction V n2mjn > 2n ' now consists of 
two terms: 

V n2n3 , nWs (p, p',E) = V%l tn , „, (p, p',E) + V^ n , n , (p, p', E) (24) 
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where we add to the potential in Eq.(p} another kernel, of the following form (cf. Eq.(26) of [l7|). 

V (QQ) _./ p\ _ ST J nn 2 n 3 (P) J nn' 2 n' 3 (p r ) 

V n 2 n 3 ,n' 2 n^P ^> - 2^ E - M n ' ( 5) 

n 

Here J(p) is obtained from J^(p,k) putting k = and replacing by M^yi23> where j/123 is given in Ea for 
different states, and M w is a fixed parameter for all charmonia and bottomonia states used in 17|-[26j], M w = 0.8 GeV, 
in actual calculations one reduces the fully relativistic vertex j/123 to the two-component spinor form, convenient 
for nonrelativistic form of participating wave- functions M^ym = 72/123, an d 7 = 1.4, see appendix C of [l7j for 
details. 

The resulting equation for rj v (E) has the same form as in (|13l) . but now A„ and K n are columns and Q is a matrix 
in n and n! indices. Fixing n' and denoting a single channel 712^3 = n' 2 n' 3 = 1, one arrives at the equation 



C12 



where 



det 1 =0, (26) 

Cl2 = C21 = J (2tt) 3 " 3 ^ Xnn 2 n 3 (p) Jn'n 2 n 3 (p) 

= E /^^(p^^nsfP)- ( 27 ) 

Solution of Eg. ([2"B]) for ^(-E = = 1 gives the pole positions E p . In particular, one can derive how the original 
(QQ) P ^ is shifted due to two effects: 1) CC to the sector II of two heavy-light states (Qq) (Qq) 2) CC due to the 
hadroquarkonium states - sector I. 

It is clear in (|2l>| , that the resulting r\ v (E) contains threshold singularities from sectors II and I and the pole at M n 
in the limit of small CC. In the spirit of calculations in [l8[ and using (f2"2"j) one can define the probability of transition 
from sector I to sector II as the absorptive part of t (|2"21 due to sector II. It yields 

\l-r) v (E)\* 

The total width T of the resonance, originating from the state (QQ) n ' in sector III is obtained from the expansion 
(l20l) for small T, and from the position of the pole E p in the equation f], y (E p ) = 1, where r\ v is found from (|26[) . for 
arbitrary 7. The partial widths of the resonance, originating from sector III, corresponding to channels in sectors I 
and II are proportional to absorptive parts of r\ v (E) on the cuts, starting from thresholds of sectors I and II. The 
concrete examples of X(3872) connected to ujJ/ip and pJ/tp states will be given elsewhere. 

IV. PION-QUARKONIUM RESONANCES 

As a specific example of the general formalism in section 2, we consider here the pion interaction with heavy 
quarkonium. To simplify matter we shall not use WEM in this section, writing all expressions in standard form, since 
we shall not use the notion of resonance wave function for pion quarkonium. The corresponding interaction term 
is given in (1), while Jnn 2 n 3 is defined in (2). However now, in contrast to the vector hadron case 
of ©, one has instead for pion emission the same vertex, which was derived in [UIH, 25]. For DD* or BB* one has 

* = X7r (29) 

while for D*D*(B*B*) one obtains 

Note, that indices i, k, I refer to the polarization states of initial and two final vector mesons. Also, the shorthand 
notation DD* implies -± (DD* ± D*D)i for isospin /state. Finally, as in ©, ^ the factor ^ is taken into 

account in the pion phase space integral over A^h . 
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With these s-wave-type kernels (l29l) . (|30[) and SHO wavefunctions in ([2]) one can write the factorized form for 

71712773 

(P- k ) = ^ ^l, ( k ) X»n 2 „ 3 (P) , (31) 

where 

vS&bOO = e"^, x™ 2 n 3 (p) - y^e"^- W 3 (cp) (32) 

and I n n 2 n 3 (p) for SHO functions for heavy-light mesons with n 2 = = 1 (which gives 95% accuracy for B,B* and 
D,D* mesons (Hll Hj|) is 

'-<■» - ^ •(-(-!). f. <»> 

where all constants are defined via SHO parameters of wavefunctions, participating in the overlap integral @, see 
Appendix 2. 

We note, that <y9^„ 3 (k) = ip(k) with a good accuracy does not depend on n 2 n 3 , when n 2 , n 3 run over a pair of 
indices B, B* (or £>, £)* etc.), since the corresponding wavefunctions are very similar. Therefore the kernel K n (E) in 
(|18|l also does not depend on indices n 2 n.3 and can be written as 

k 2 

/<i 3 k e 
(2tt)3 2^(k) (S„(k)+o; ff (k)-£;) ' (34) 

Defining as in (|14p 

> / p \ _ _j_ /" ^ 3 P Xnn 2 n 3 (p) Xri'ri 2 n 3 (p) /or\ 

^ ~ V (2-) 3 E n2 (p) + ^ 3 (p) - e ■ 

one has a system of equations (|17p. from which one defines resonance energy [l7j ] 



det 



l-*f(£)C(£) =0, (36) 



where (#)„„< = K n 8 nn > and (CW = Cnn'- 

Note, that Xnn 2 n 3 depends on polarization states of all particles, and that of n (index i in (|29p. (|3"0)l ) can be fixed, 
while one should sum up in (|35j) over spin and isospin projection of particles ri2, ns- 

At the end of this section we consider the contribution of pion-quarkonium resonance into production crossections 
of the final state (QQ) n ' ^it. In the zeroth approximation the amplitude for the transition {QQ) n {QQ)ri tttt was 
calculated in [24l426l |. 

„,(iir)/m \ " d 3 p Jnn 2 n 3 (p,ki) J m „ 2 n 3 (p,k 2 ) , - 
k 



E 



(2tt) 3 


E 


-£n 2 n 3 (p) -^(ki) 


d 3 p 


J (2) , 

nn 2 


rl ,( P ,k!,k 2 ) J^ n , n ,( P ) 


(2tt) 3 


E - 


£„ in -(p)- J E(k 1) k 2 ) 


d 3 p 


J nn'2 


<(P) J m<n 3 '(P' k l' k 2) 



fi/tg ''■3 V"/ 771 71 2 1^3 * ' (0*~7\ 



Here Jnn 2 n^ — ^77712713 (p; k) and 



(k 1+ k 2 ) 

J^Cp^k^-^^e ^ e-^r/ TO2R3 (cp); (38) 



8 








FIG. 3: Rescattering series yielding a possible pole structure in (T(n')7r). (Contribution to a nm ) 




FIG. 4: First order (in contributions to the factor b nm in Eq. (|43|l . 



with 

123 /I v /2c 7r (k 1 )2 W7r (k 2 )^ ' 

while Jrm 2 n 3 (p) is the pionless overlap integral ([2]), where k = 0, and ~~ * Z/123 = M u — ((ft — 2 f P '+ 9 — 7)- 
Looking at (|37|l . one can realize, that it can be written as 

(2) w ( ™\E) (40) 

and the first two terms of ^w^m\E) are 

My,™ = ^(kx) C„ m (^') p(k 2 ) + p(k 2 ) ( nm (E") ^(kx) , (41) 

where ( nm depends on the energy E' and E" of the (QQ) 7r 2 and (QQ) tti systems respectively, its Lorenz invariant 
definition see below. It is clear that these terms are the first terms of the whole rescattering series, depicted in Fig{3j 
which can be summed up as follows 



(D w (-) 



nm 



-+ (1) w nm = ^k 1 ) (c— y <p(^») + ^2)(cyzkc) ^ k i)- ( 42 ) 



Note, that for the transition T(n) — > Y(n')inr, C(2S') and K(E') in the first term in (|42)) depend on the invariant 
mass of T(n')7r 2 , while in the last term on the r.h.s. of (|42l) . Q{E") and K(E") depend on the invariant mass 

Mfol of T(n') 7Ti, 

M ml = ^ ~ 2 A/nWi + mj, i = 1,2 

We now turn to the last two terms in (|37j) . which contain two-pion vertex, shown in Fig|4j and take into account, 
that the dependence on ki,k 2 there is contained in the factor <p(ki + k 2 ), as was shown in [24|, [2f|, as well as in 
the series shown in FigJS] (and the similar one with ki f-> k 2 ). Hence the total amplitude of (n,ni) transition with 
emission of two pions can be written as 

Wnm = M nm = <p(ki) <£(k 2 ) a nm - (p (ki + k 2 ) b nm (43) 
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FIG. 5: Rescattering series including double pion production vertex. (Contributing to 6„m)- 



where 



C(£') 
C(^) 

One can see, that (|43|) satisfies Adler condition 



1 



1 - K(E') C(E') 
1 



(£' -> E") 



jr + (E'^ E"\ 



1 - K{E')C(E') 



a nm {h = 0) = &„ ro (A;j = 0), .M nm (ki = 0) = 0, 
and for not very large ki , k2 one can write approximately 

M (^(ki)p(k 2 ) -y?(ki +k a )). 

The probability of transition (QQ) n — > (QQ) n ' 7r7r is 

dw((ra) -> (n'jTTTr) = |7W„„'| 2 f^a ^fr^' +^(ki) + ^(kj) E„ ) 



(2tt) 3 (2tt) 3 2wiw 2 



and the dipion decay width is 



t(7T7t) 



dw((n) -> (n')7T7r) = J d§ \M r 



where is the phase space factor 

1 (M 2 + M 2 , - q 2 )(M n + M n ,) 



32 TT 3 



4M 3 



- V( A -^) 2 - 9 2 Ve 2 - 4 ™ 2 dq d COS V 



with the notations 



AM = M n - M n > 



g 2 = M 2 OT = (h + k 2 f = {w t + lo 2 ) 2 - (kr + k 2 ) 2 . 
Finally, one can also study the process e + eT — > (QQ) n ' tttt, with the amplitude 

1 



E + w-E 



M nn > 



(44) 

(45) 
(46) 

(47) 
(48) 
(49) 



(50) 



(51) 



where M. n n' is given in (|43| and c m = 47rc y 2 Qv/ ^ -0 m (O) so that the contribution AR of (QQ) n > tttt to the hadronic 
ratio R is 



72 7re^ 

Afl((n')7nr) = -^ 



Y>n(0)U " ) M mn ' 



(52) 
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TABLE I: List of thresholds (in GeV). 



Threshold E th Threshold E t 



h 



BB* 10.605 T(lS)n 9.60 

B* B* 10.650 T(2S)tt 10.16 

B S B* 10.780 T(3S)tt 10.495 

B*B* 10.830 T(4S)tt 10.720 



T(5S)tt 11.00 




FIG. 6: Real (solid line) and imaginary (dashed line) parts of (33(E) and K3(E) (see Eqs. (|34ll35p ). One can see cusp structures 
at the thresholds BB*, B* B* and Y(3S)tt correspondingly. 



We now turn to an example of possible pionic bottomonium state in the reaction T(55) — > (T(n') tt) tt, which can 
proceed through the chains 

T(5S) -> J2 ( B BUn 3 n -► ^(T(n') tt) tt. (53) 

riin^ n' 

We are interested in possible poles in the J p = 1 + channel of connected (BB) n2n3 and T(n') 7r states, which are given 
by the equation 

det [<W ~ K n {E) C nn > (£)] = 0. (54) 
Neglecting first nondiagonal elements of Cnn' ; one has an equation for E 

K rmc (F)- 1 V f d3p XnmnM f d 3 k ^(k) 

" l J JVcf^y (2tt) 3 S n2 (p) + ^ 3 (p)-^7 (2 7 r) 3 2 W7r (fc)(K(k)+u; 7r (k)-£) 



One can easily recognize in ([55]) the norm of the kernel of the integral equation . 

The analysis of (155]) starts with the list of thresholds in (BB) n2Tl3 and in T(n') 7r channels in Table fl] One can see, 
that the most important combination is T(35) tt O £?£?*, B*B* with additional contribution of T(25) tt, T(1S) tt and 
B S B*, B*B* as a next step. Since the maximal energy of our systems in the reaction ([53)1 is M(5S) — — 10.725 
GeV, one is mostly interested in the energy region 9.78 GeV < E < 10.72 GeV. 



V. RESULTS FOR T(5S) T(n'S) nn 



In this section we study numerical results for the reaction T(55) — > T(n'S) tttt, and we shall be interested mainly 
in the possible appearance of resonance-like structures in the systems T(n'S) tt, n' — 1, 2, 3. The resulting equations 
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q\ , -^j , 1 o 1 — — ' 1 1 1 1 — 1 ' 1 

10.5 10.55 10.6 10.65 10.7 10.5 10.55 10.6 10.65 10.7 10.5 10.55 10.6 10.65 10.7 

M(Y(1S)ti) [GeVj M(Y(2S)7t) [GeVj M(Y(3S)tt) [GeV] 



FIG. 7: Modulus squared of the first term in (|44|l |a nn /| 2 as a function of the invariant mass of (Y(n'S)Tr), computed for the 
reaction T(5S) -> T(n'S) tytt, n' = 1, 2, 3. 

for differential and total probabilities are given in Eqs. ((33]), (fl6")l , (H71) . The coefficients y[^\ and parameters of 

T(nS) and £?, £?* wave functions, needed for calculation of £ nm are given in Eqs. (|29l [30]) and Appendices 1 and 2. 

It was assumed above, that the knowledge of wave functions and channel coupling constant M u (one for all types 
of strong decays) can describe all CC phenomena and, in particular, level shifts due to the CC. However, at this 
point one encounters the fundamental difficulty, which was studied in [36| and is not still resolved. The point is, that 
assuming a constant CC, not depending on participants of decay process, one obtains huge shifts of energy levels 
(several hundreds of MeV) due to mixing with higher states. (To improve situation, in [l8| a cut-off coefficient k ps 0.5 

was introduced for contribution of higher levels). This calls for a detailed scrutiny of our matrix elements Jmon,(P : k) 
and basic matrix elements (without hadron emission) Jnn 2 n 3 {p), entering in the expressions for energy shifts, see [1 71 ] 
for details. Indeed in Eq. ([2]) one can see, that the string in the original hadron n, placed between points u, v is 
decaying at point x into two hadrons, placed between u, x , and x, v respectively. It is clear, that when the point x is 
far away from the center of the string u, v, the string breaking process does not occur. Replacing distances between 
points by typical radii i? n , R n2 > Rn 3 of states n, rt^ one obtains condition (for R n2 sa i?„ 3 ) i?^ 2 < jR^ + p 2 , where 
p is the typical string width, p w 2A w 0.3 fm [371 1381]. The corresponding factor can be rigorously deduced in the 
formalism of (28|, and the resulting cut-off strongly decreases the CC between states with incomparable sizes. This is 
taken into account below by assuming different M u for different decay channels. We take for M u the value 1.0 GeV in 
case of n — 4, 5. For n= 1,2,3 we take = 0.1, 0.2 and 0.3 GeV respectively, to take into account decay mismatch 
between the sizes of T(nS) and BB* systems, since R(IS) = 0.2 fm and R(2S) = 0.4 fm, while R(B) m R(B*) ps 0.5 
fm. This mismatch is not taken into account for simplicity reasons in the general definition of the overlap integral 
©■ 

In the beginning we have estimated K n (E) and ( n m(E) in (]55j) for ji, m = 1,2,3,4,5, using parameters of SHO 
functions, which were used before in [25L l2q|. they are given in the Appendix 2. Real and imaginary parts of C33(i?) 
and K 3 (E) are given in FigJH] The possibility of peaks e.g. in T(n'S)w system is demonstrated in FigJT] where we 
plot the quantity |a5n' | 2 for n! — 1, 2, 3 respectively (or more exactly the first term of (|44]l ) as function of the invariant 
mass M-^j. One can see sharp peaks near the thresholds BB* , B*B* at 10.6 and 10.65 GeV respectively. In the total 
distributions, however, a more complicated combination of terms enters, as seen from (|44[) , and one should calculate 
a symmetrized in M^l, M^ v expression, and moreover take into account appropriate phase space factor. One can 
associate these peaks with poles, situated in the vicinity of these thresholds (see discussion in Appendix 3). 

Note, that two peaks in FigJT] occur from combination T(n' S)iti. However in the full decay distribution du>((5) —¥ 
(n')TTTr) the symmetrized sum l|44p enters, which produces an additional mirror reflected pair of peaks in dw if 
considered as function of M^ v , q, since 

(M&l) 2 + (mJ$) 2 = Ml + Ml +2ml-q 2 . (56) 

The full probability distribution (|47l) contains symmetric sum of two rescattering series for T(n'S) 7Ti and T(n'S) 7r 2 
respectively. When plotted as function of M^ v it contains both poles of the first series at M^J V — 10.610 and 10.650 
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0' ' * ' 1 1 1 — 1 1 1 1 ' 1 ' 1 

10.5 10.55 10.6 10.65 10.7 10.5 10.55 10.6 10.65 10.7 10.5 10.55 10.6 10.65 10.7 

M(Y(1S)ti) [GeVj M(Y(2S)tt) [GeVj M(Y(3S)tt) [GeV] 



FIG. 8: The distribution dT „„/ / dMi„ v for the transition T(nS) — > T(n'S) mr as a function of the invariant mass of (T(n'S)n) 
for n' = 1,2,3. 

GeV, and also poles from the second series at the points M^ V (M^ V = 10.610, 10.650) defined by Eq. (f56|) . In this way 
one obtains Fig|51 where for the case (n, n') — (5, 1) and (5, 2) the secondary poles are out of mass interval, while for 
the case of (5,3) the pair of secondary poles overlap with the proper poles of M^ v . In experiment one can separate 
points on Dalitz plot relating to M-^ v and M^ v , which results in two plots with the same pai r of peaks. 

Recently experimental data on the reaction T(55 l ) — > T(n'S) itit, n' = 1, 2, 3 appeared in [34| . and the experimental 

distributions duo ((5) — > (n'Jmr) are presented in Fig. 5 of [34| as functions of min{M^2, M-^ v } and max{M^2, M^ v }. 
In order to compare our results with experimental data, we calculate distribution (|47p in terms of invariant masses of 
T(n'S) 7T and tttt systems. We use expressions (|46p . (|4TJ)) and formulas from [24[ to express quantities like w ff (ki), (ki + 
k2) etc. via variables (q, Mi nv ). After integration over q we obtain the distribution in terms of Mi nv (T(n' S) ir). The 
result is presented on FigJS] Comparing experimental data (left and right plots, glued together for n' = 3 and left 
plots for n' — 1,2 on Fig. 5 from [3J]) with our theoretical calculations, one can see a close similarity in the general 
form and position of Zf,(10610) and Zf,(10650) peaks, which in both cases appear at BB*, B*B* thresholds for all 
available n' — 1,2,3. We do not intend to reproduce here all features of experimental data, which depend strongly 
on details of wave function profile, but we study mostly the dynamics of process. For the better description we only 
slightly changed w.f. of T(nS), presented in Appendix 2, by ~ 20%. A more detailed quantitative comparison of our 
data with experiment (34|, and decay distributions as functions of q = and cos#, are now in progress and will 
be reported elsewhere. 



VI. CONCLUSIONS AND OUTLOOK 



Summarizing our results, we shall stress the main features of our approach. We have considered the interaction 
of a light hadron with heavy quarkonium, arising solely due to transition to intermediate states of two heavy-light 
mesons. No direct interaction between light hadron and heavy quarkonium or else between two heavy-light mesons is 
assumed, therefore our dynamics has nothing to do with molecular states in the strict sense. Our mechanism is also 
distinct from dynamics of (4g) or hybrid states. 

We have constructed explicitly transition vertex for the strong decay with emission of a light hadron from the 
first principles, and then all dynamics is defined by the overlap integrals of all wave functions involved, i.e. heavy 
quarkonia and heavy-light mesons. The latter have been found in 35] from the relativistic Hamiltonian, containing 
only first principle input, and accuracy of Gaussian representation was checked in [i~7l l32l|. 

The answer given in Figs l7l8l is positive for the case of pion-bottomonium system, and agrees with the recent exper- 
imental data [34| at least qualitatively. More detailed calculations are still needed to check all results quantitatively 
vs experiment and this program is now under study. Other systems should be treated as well, e.g. nhb, studied 
in [3J|, and the series of Z(4430), Zi(4050), ^(4250) resonances in the (cc)tt system. All formalism used above for 
bottomonium can be applied to the charmonium case without modifications; for the vector charmonium or bottomo- 
nium (/if,, Xb, h c , Xc) one can use transition vertices, given in Appendices 1,2, and symmetry properties of the whole 
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amplitude will be different. This work is now in progress. 

One should stress at this point, that resonances, found above in T(n / S)tt system are specifically multichannel ones, 
in the sense, that they belong (and appear) in all n' channels and are given by zeros of det[l — KQ. 

At this point it is convenient to compare predictions of molecular and CC models for the positions of resonances. 
As one can see in Figf?! our CC model predicts sharp peaks at BB* , B* B* thresholds (possibly due to virtual 
multichannel poles) in all Y(n'S)ir (n' — 1, 2, 3) channels, and this effect is due to combined contribution of £ and K 
terms, i.e. pure CC effect. In contrast to that, the pure molecular picture, when poles in the BB*, B*B* systems 
appear due to direct internal interaction (i.e. without irT(nS) channels), the poles appear in all £ nn (E) terms near 
corresponding thresholds. This can be seen simply in the definition of C(E) in (|35p . where XnXn' is multiplied by 
the free BB* (B*B*) Green's function. In case of strong interaction in BB* , B*B* systems, this Green function is 
replaced by the exact one and should contain pure molecular poles i.e. C, nm (E) = . Insertion of this form into 

the det[l — KQ = yields the n-th order equation for roots in energy, and these n roots are all strongly displaced 
from the original places at the thresholds for large K( coupling (large C n ), and are almost degenerate n poles at 
thresholds for small coupling. Both pictures are different from the experimental data - the same two poles at BB* 
and B*B* in all irY(nS) and irh^mP) channels. Thus the situation with the only pole at each threshold is possible 
and characteristic for our multichannel CC resonance and is unlikely for pure molecular states. 

Another property is, that resonances appear most likely, when thresholds in K n and Q nn for some n are close to 
each other, and then this channel n will be the dominant one, making det[l — KQ close to zero. In the case of 
%Y(n'S) the dominant channels are those with n' = 3, 4, as can be seen from Table in section 4. Therefore the visible 
channel h (QQ) n >, where a peak is found, is not necessarily the dominant one, as might be in the case of Z(4430) with 
a peak seen in TT + ip' channel. The dynamical reason, why proximity of thresholds is favorable for the appearance of 
a resonance, is that both Ke[K n (E)] and Re arc decreasing fast enough away from thresholds, making their 

product maximal for the coinciding thresholds. 

In all discussion above the notion of resonances (or virtual and real poles) was stressed, and hence the whole sum 
of rescattering series as in Figs l3|5[ were considered. But it is possible, that already the first terms of this series can 
contribute to enhanced correlations, which look like bumps in decay distributions. This approach was considered in 
two recent papers [H|,[4(|, and is especially important in the case, when high spins and angular momenta are involved. 
Relation of this approach to our methods of present paper seems to be practically important and the corresponding 
work is planned for the future. 

Recently several papers appeared U 42|, where i?b(10610) and Zf,(10650) are considered as molecular states and 
treated in the QCD sum rule method [4l|, while in [42| these states were originally supported to be xh and Xb ± 
shifted to the BB* and B* B* thresholds (as it happens in A (3872) case). As was argued above in both cases the 
poles produced appear in £ nn ' and would be strongly shifted in the rescattering series of FigJU yielding 5 peaks at 
different masses in Y(nS) n, hb{mP) tt, n = 1, 2, 3; m — 1, 2. 

In [43j an interesting analysis of T(55) — > Y(2S)irir decay is presented, demonstrating important contribution of 
Zf,(10610), Zb(10650) to the decay distribution in M(tt + it~) and cos# a detailed comparison of these results with our 
approach and previous results in [26] is now in progress. 

The authors are grateful to A.M.Badalian, R.Mizuk and P.N.Pakhlov for many useful discussions. The financial 
support of Dynasty Foundation to V.D.O. and RFBR grant no. 09-02-00 620a is gratefully acknowledged. 



Appendix 1 
Calculation of the transition kernel y[*£ 3 



We use here, as well as in [17|, [24], [25[ the spin-tensor representation of wave functions and Z factors instead of 
technic of Clebsch-Gordan coefficients. We start with the (4 x 4) fully relativistic technic given in @, H^], where §123 
can be written through the so-called Z factors (we omit the superscript h for time being, cf. Eq.(A2.24) of [24]) 

Z%23x , , , , x 

2/123 - 7= (Al.l) 

where Zyi^x and Zi are Z factors for the total process and for individual hadrons respectively participating in the 
transition process 1 — > 23 or 1 — > 23 h, shown in FigJU Defining projection factors for quarks and antiquarks Ajf, 
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TABLE II: Bilinear operators tpTitp and their (2 x 2) forms (Notations see in the text). 



(2 x 2) form 



o-+ 

\— 
1+- 
0++ 
1++ 
2++ 
2- + 
2 

1"" 



3 5r 
'ft 



3 Pi 

3 £ 2 



-«75 
Ji 

-475 <9* 

1 

7i <9fc +7fc 9i -f &fc# 

-|5 lfc (9) 2 )75 

(7« dk +7fc 9i -f *ifc9)7B 

7i^ife 



1 

v-2 



0» 

crn 



V V — w w 

— (v c aiV + w c aiw) 
v c di v — w c di w 

i(v c w — w c v) 

— (v C <JiW + W C CJiV) 

— (v c pikv + w c p ik w) 
i(v c Lo ik w — w c uj ik v) 
~(v c pikW + w c p ik v) 
— (v c aiUJikV + ■w c a i uji k w) \oi(njn k — ffa) 



75"- 
^■e lk ia k ni 



Here P2(o", n) = §(ff>raj + a e rn 



(an) 5 U ). 





3 

FIG. 9: Decay matrix vertices. 



where subscript fc refers to light quarks, k = q, or heavy quarks, k = Q, one can write 

^123* = tr(r x a+ r 2 a- r x a+ r 3 Aq), 



(A1.2) 



= tr(Ti A+ r, A" 



and 



i ■ (fe) 

± _ m k ± u kli T 7* , _ „ n 



Note, that the sum over (i) in (A1.4) is for i = 1,2,3, also in the cm. system p^' = —pY' = Pi, while m k is the 
current quark mass and uj k is the average kinetic energy of quark k in the hadron. 

The operators Ti correspond to quantum numbers of a hadron, and are given in Table 1 below, while T x refers to 
the process under investigation, for the case when no hadron is emitted, T x = 1, while for the pion emission T x ^ = 75 
and for vector particles T 
In this way one obtains 



,(2) 



(A1.3) 



(A1.4) 



(«) 



7»- 



J D,B 



n 



4 02 
3 S V.T 



■ mz 



n 2 



(A1.5) 



where £Id b 



"c,b, 



D,B, 



n 



ni 2 h )v,x are calculated in [24j and given in Appendix 1 of [14]. 



Examples of relativistic (4 x 4) expressions of 2/123 are given in [2J, |25|, e.g. for (1" 
(1")„ -> B*B*n (D*D*tt), y[l\ are given in p9 ]) .(|30 | while for (1~) -> BB* 



BB*n (DD*tt) and 
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We now turn to the case of (2 x 2) formalism, introduced in pjj, where resulting kernels are denoted as 2/123, and 
are computed according to 

vlt - *r{r&> rt=> ( ffq ) r£*>} (Ai.7) 

while for emission of an additional pion in y^l one should omit the factor (erq) in (|Al.7|) . In this way one 
obtains (|29|) . ([30]) . Note, that normalization of states in (2 x 2) formalism is different, and one should sum up 

over all polarizations in initial and final states (extra factor of l/\/2 is in (IA1.7|) as compared to (|A1.6I) . and — — > 
1 " 

m+e n + {U)-(V) ■ 

Of special interest are the transition kernels for P states, where one has for 3 Pi — > (BB*ir), y^l r ° d = ^ Zini ni 
and for 1 P 1 -> (BB*ir), y[f 3 rcd = 0, while the (4 x 4) formalism yields 

2/123 = —T TT ~ eklm ( k ~ ( 1 / 2 P~ 1)1 ■ ( A1 - 8 ) 

4 ilQ p 

For 1 P 1 -> (B*B*ir) one has rod - n t 6 kl , and for 3 i\ -> (5*5*71), j/W rcd = » ^ n, - 5 si 5 ti ) i, s, t 

are indices of polarizations. 



Appendix 2 

Wavefunctions of heavy quarkonia and heavy-light mesons 

In Eq.© R<gf and are series of oscillator wave functions, which are fitted to realistic wave functions. 

We obtain them from the solution of the Relativistic String Hamiltonian, described in [35| . 
In position space the basic SHO radial wave function is given by 

R^W, r) = ^ T{ l+i + y2 )^y e_/3V/2 ^-/W) (A2.1) 
(R s n f°(P, r))Vdr=l 



1 



where f3 is the SHO wave function parameter, and L 1 ^}^ ((3 2 r 2 ) is an associated Laguerre polynomial. The realistic 
radial wave function can be represented as an expansion in the full set of oscillator radial functions: 

R nl {r)=Y.c k Rtr{P,r). (A2.2) 

fe=l 

Effective values of oscillator parameters (3 and coefficients c k are obtained minimizing \ 2 and listed in the Table HIT? . 
In the momentum space the SHO radial wave function is given by 3 



H m [P,P)~ ^3/2 Vr(n + / + l/2) \j3 ' 1 m 







(i^ r ( / J,p)) 2 g*=l (A2.3) 



2 Typos of the sign convection are corrected for the IS, 25, 35 charmonium states of |17| 

3 Note a typo in the equation for B^°{fi, p) of [l]} 
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TABLE III: Effective values /3 (in GeV) and coefficients Ck of the series of oscillator radial wave functions R^." (P, r) which 
are fitted to realistic radial wave functions R„i(r) of charmonium, bottomonium, B and D mesons. 



State 


P 


Cl 


C2 C3 


C4 


C5 








Bottomonium 






IS 


1.27 


0.977164 


-0.151779 0.141319 


-0.020857 


0.036863 


2S 


0.88 


-0.18823 


0.953901 -0.135824 


0.181277 


-0.001509 


3S 


0.76 


-0.128081 


-0.145885 0.936255 


-0.169962 


0.226281 


4S 


0.64 


-0.019432 


-0.149876 -0.362548 


0.88647 


0.014308 


5S 


0.6 


-0.011183 


-0.016911 -0.182019 


-0.403138 


0.853936 


IP 


0.93 


0.977994 


-0.165514 0.122975 


-0.018631 


0.024364 


2P 


0.76 


-0.092083 


0.971982 -0.137347 


0.162174 


-0.002707 


ID 


0.8 


0.979042 


-0.168619 0.11132 


-0.016998 


0.018383 


2D 


0.69 


-0.063135 


0.979871 -0.117356 
Charmonium 


0.145358 


0.000973 


IS 


0.7 


0.97796 


-0.169169 0.117682 


-0.019694 


0.025113 


2S 


0.53 


-0.121144 


0.973054 -0.130808 


0.141495 


0.00097 


3S 


0.48 


-0.096897 


-0.086156 0.961504 


-0.155931 


0.178935 


4S 


0.43 


-0.021639 


-0.128342 -0.215657 


0.947215 


-0.028308 


5S 


0.41 


-0.010701 


-0.022826 -0.157891 


-0.258875 


0.925921 


IP 


0.57 


0.976869 


-0.184163 0.105506 


-0.018941 


0.017215 


2P 


0.48 


-0.063059 


0.981868 -0.123012 


0.127035 


0.000588 


ID 


0.51 


0.979118 


-0.178313 0.095356 


-0.016095 


0.013123 


2D 


0.45 


-0.044316 


0.986084 -0.107408 
D meson 


0.117002 


0.001907 


IS 


0.48 




c=l 
B meson 






IS 


0.49 




c=l 







Then using Eq. (A16) from Appendix A of [25(, one can write I n n 2 n 3 (p) m (E2) from = = 1 as (/3j = 

ft»,A = A») 



where 

A n = 2/3 2 1+ f3l y= Zl , ?2 > f = -t-^-^t=j (A2.5) 



2/31 - 0j 2pft 



2vr\ 3/2 (2^//3 2 ) 3 



Here $ is the confluent hypergeometric series, 



Pi J v /2 2 » 7 r 3 / 2 (2n- 1)! 



a a(a + 1) 2 



$ (a;7 , x) = l + _, + ___ x 



(A2.> 



Appendix 3 

Analytic structure of hadron-qurkonium resonances 
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Re[K(E) 




\Ethr 



FIG. 10: a) Contours C and C' in the complex E plane, exibiting analytic properties of the integral (IA3.4[I : b) Typical behaviour 
of the Re[K n (E)] and lm[K n {E)} given by Eq. (MM- 



At the end we study the analytic structure of the resonance denominator in - ^ ~ , which in the one-channel case 
is given by Eq. (|55|). One can write K n (E) as the integral 

W^f^, (A3.) 

where k(uj) = yjij 1 — m\. To display the analytic properties of K n (E), we shall use the method, exploited in pj}, 
Appendix E, i.e. we replace the integral in (|A3.1I) as 1/2 of the contour integral over the contour C, circumjacent to 
the interval [rn w , 00], as shown in FigflUl and take E in the physical region above the contour C, E — s- E + iS, 5 > 0. 
Deforming the contour, so that it passes above and to the left of the point E (contour C in FiglTUl). using Cauchy's 
theorem, one has the representation 

K n (E) = ± <p 2 (k(E)) k(E) + F K {E) = i a K {E) k{E) + b K (E). (A3.2) 

47T 

Here k(E) = \J{E — M n ) 2 — m%, and Fk(E) is analytic function in the neighborhood of the threshold — 
M n +171^. Hence Re[K n (E)) acquires a negative contribution from the first term on the r.h.s. of (|A3.2I) and behaves, 
as shown in Fig[inj This behavior of Ke[K n (E)] and lm[K n (E)] agrees with that, obtained by numerical integration 
in Fig|6l 

In a similar way one can write the form of £ nm (E) 

( nm (E) = £ 4 P j^ e / jf (P) \ = i^T^W x'W*)) + F C^} =iac{E)p(E) + b ( (E) (A3.3) 

\2M ~ Z 



where p(z) — v2Mz, z = E — M„ 2 — M„ 3 , M = m™ 2 +m 3 ' anc ^ ^ (* = C) are analytic and positive functions 
near the thresholds. 

From (|A3.2[) , (|A3.3[) one can deduce, that the product K n (E) ( nn (E) , is a real analytic function in the complex 

(n) 

plane of E with cuts, starting at thresholds E^' — M n + m^ and E n2 n 3 = Af„ 2 +M„ 3 and going to plus infinity. Thus 
the combination (1 — K n {E) Q nn (E)) is a real analytic function in the E plane with positive weights in the integrals 
(IA3.2|) . (|A3.3I) . Hence the only possibility for the zeros of (1 — K n (E) ( nn (E)) is on the real axis below thresholds, or 
else on the next sheets, which implies a standard situation with possibility of bound state or virtual state poles, or 
else Breit-Wigner poles in 1 _ K ^4-^ — t^j. From (IA3.2|) . (|A3.3p one has 

(A3.4) 



1 - K n (E) (nn(E) 1 - b K fe c + kpa K olq - i (ka K b^+pa^ b K ) ' 
In the simple case, when the thresholds in K n (E) and Cnn(E) coincide the resonance factor acquires the form 

1 _ 1 

1 - K n (E) Cnn(E) ~ A-ikB 



(A3.5) 
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with B > 0. This form demonstrates the appearance of a virtual (A > 0) or a real (A < 0) pole. 
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